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Abstract 

The two-dimensional random-bond Q-state Potts model is studied for Q near 2 via the 
perturbative renormalisation group to one loop. It is shown that weak disorder induces 
cross-correlations between the quenched- averages of moments of the two-point spin/spin and 
energy/energy correlation functions, which should be observable numerically in specific linear 
combinations of various quenched correlation functions. The random-bond Ising model in 
(2 + e) dimensions is similarly treated. As a byproduct, a simple method for deriving the 
scaling dimensions of all moments of the local energy operator is presented. 

1 Introduction 

When studying disordered systems at or near their critical points, the notion of a single scahng 
dimension governing the behaviour of an observable field at large distances becomes inadequate. 
Randomness gives rise to a broad probability distribution for the scaling dimension of a two-point 
correlation hmction InG/ Inr, so the quenched average of n"^ moments of the two-point function 
will have a non-linear dependence on n: 



G"(0,r) r"^^" with x„ ^ nxi 

where the over-bar indicates the average over quenched disorder. Such multiscaling behaviour has 
been both predicted analytically and observed numerically in a wide range of systems, but we will 
primarily concern ourselves with the random-bond Q'State Potts model in 2D (the random-bond 



Ising model in (2 -|- e) dimensions is briefly considered in section 3.4). For Q near 2, the critical 
behaviour of this system is obtainable via an e-expansion about the pure Ising- model fixed point, 
with an expansion parameter proportional to the specific heat exponent a. (We recall the Harris 
criterion |l[], which states that bond randomness will be relevant if and only if the specific heat 
exponent is positive: hence the disorder is marginal for the Ising model ((3=2) and relevant for 
all Q greater than this value). Following initial work by Ludwig [^,^, multiscaling behaviour of 
moments of the two-point spin-spin correlation function has been well established in this model 
both theoretically and numerically (see e.g., |4|-|^), and multiscaling in the corrections to scaling 
for the energy-energy correlator has been predicted. (The energy-sector calculation is complicated 
by the fact that scaling operators correspond to irreducible representations of the permutation 
group of replica indices, Sn)- 

In this paper we examine some further structure induced by the disorder, which shows up 
as correlations between the two-point spin/spin and energy/energy functions. (In fact, the non- 
trivial replica structure of the energy sector means that a distinction must be made between the 
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connected and disconnected two-point energy correlators: these may have different behaviour with 
the spin/spin function, and also be mutually correlated). For example, does 



(a(0)a(i?))f(£(0)e(i?))9(,)^i?- 



or is a more complex behaviour to be expected? We show that this is dependent on both the 
model in question, and whether the connected or full (ee) correlator is taken. 

In section ^ we shall introduce the Q-state Potts model and the replica method used to carry 
out the average over quenched disorder, before describing results already obtained for this model. 
In section |3.1| we present calculations of scaling dimensions (to one loop) of moments of the 
energy-operator, and then show (section [3.2| )that a non-trivial behaviour of the mixed spin-energy 
moments is to be expected. This indicates that the quantities In Gaa/ In R and In Gee/ In R should 
be considered as being drawn from a joint probability distribution function that is not equal to the 
product of the individual marginal distributions, i.e., the two quantities are not statistically inde- 
pendent. We then make contact with quantities available via numerical calculations, by expressing 
two-point functions of the irrep scaling operators as linear combinations of various quenched av- 



eraged energy correlation functions (section 3.3). Finally we describe a similar calculation for the 



random-bond Ising model in (2 + e) dimensions (section B.4) 



2 Model and previous results 

We write the reduced lattice hamiltonian for the Q-state Potts model with weak bond disorder as 

fcT 

where e {1 . . . Q} are the Potts spin variables, Jq is the average bond strength and SJi j <C Jq 
are the local fluctuations about Jq, assumed to be completely uncorrelated. The sum is taken over 
all nearest-neighbour pairs. Taking the local energy density 

e(r) = e(ij) = ((-Jo + 5J)/kTc)Ssi,s, 

and using replicas to average over the disorder, we end up with an effective hamiltonian in terms 
of n replicas of the system coupled together. We assume that n is large, and take the n — > 
limit to perform the quenched average at the end of the calculation (see or chapter 8] for 
an overview of the replica method in field-theoretic calculations): 

n / \ n 

H,f! = Y,iH:{Q) + tY,ea{r)\-Y.Aj2ea{r)et (r) 

a=l \ r / r a,b 

where H*(Q) is the reduced hamiltonian for the critical pure Potts model, t is the reduced tem- 
perature (r — Tc) /Tc and A is proportional to the second cumulant of the distribution for disorder 
(power-counting arguments show that higher cumulants are irrelevant at the pure fixed point for 
Q near 2). We note further that the expectation value of e only shifts the critical temperature and 
cumulants, without affecting the critical exponents in any way, so we can write e(r) = (e) + e{r) 
and absorb the first term into the definitions of the other parameters. Also, terms with a = 6 in 
the double sum can also be neglected since they are either irrelevant by power-counting or yield 
disconnected diagrams which do not contribute to the renormalisation appendix B]. Finally, 
moving to the continuum limit, we find 



iJeff=^ 



a—1 ^ a—1 ^ a^b 
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Correlation functions calculated against this effective, 'replicated' hamiltonian will correspond to 
correlators averaged against the initial hamiltonian with quenched disorder: 



(e(r))ff < — > lim(£i(r)}i.op 
{e{0)eir))H ^ lim (£,(0)£,(r)},,p 

n—>-0 

{eiO))H{eir))H ^ lim (£,(0)e,(r))„p, z ^ j 

where i and j label the replicas the operators are lying in. Note that it is assumed that the above 
correlators are independent of which replicas are actually taken — this is the assumption of replica 
symmetry, which appears to be valid for weakly disordered ferromagnets. 

The quenched average of p^^ moments of the spin operator was calculated by Ludwig for all p, 
to one loop. The scaling dimension of the operator cti(T2 . . . Cp was found to be 

Xp = px, - ^p{p - 1) + Oiy^) 

where Xa- is the scaling dimension of the spin operator at the pure fixed point and y is the 
RG eigenvalue of A, vanishing proportional to (Q — 2). For typical realisations of disorder, the 
behaviour of the two-point function at large r is governed by a multifractal exponent (c.f. [^). 
This is given by the saddle point of the Lcgcndrc transformation of Xp, and was found to be 

= x, + y/16 + 0(2/2). 

In general a given moment of the energy operators does not have a pure scaling behaviour. 
In the replica formalism, our perturbing operator A X^a^^fc ^a^b is a singlet under the group of 
permutations of the replica indices, Sn (in this context the permutation group is often named the 
'replica' group instead) . Thus scaling dimensions of operators will only necessarily be constant on 
subspaces of operators corresponding to irreducible representations (irreps) of Sn- Degeneracies 
between irreps may arise — in fact all irreps of moments of the spin operator are degenerate, 
due to their different fusion rules with the perturbing operator — but in general a given operator 
6162 ■ ■ - Eq will have a 'sum-of-powers' scaling behaviour 

(£162 • ■ • eq{Q) £i£2 . . . eq[r)) - ^ ^^r^^^" 

where /i runs over the different irreps of Sn present in the decomposition of £1 £2 ... £5. 

Calculations for q < 3 were given in ||^. In each case the most antisymmetric irrep was the 
most relevant, and the scaling dimensions of these most relevant irreps were linear in 5, suggesting 
no multiscaling behaviour to this order. A more complex structure was however present in the 
corrections to scaling. 
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3 Calculation 

3.1 Energy Sector: irreducible representations of Sn and SU{2) 

First we shall exploit the connection between the irreps. of Sn and SU (2) to obtain complete sets 
of scaling dimensions for moments of the energy-energy correlation function, (£(0)£(r))'J, to one 
loop. At this order, the effect of a shift in the short-distance cutoff on the couplings is most easily 
treated in terms of the operator product expansion (OPE) (see e.g., chapter 5]). Consider a 
general partition function Z for a critical hamiltonian H* , perturbed by a set of scaling fields (pi 
with corresponding couplings 5^, 
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where a is the microscopic cutoff. The couphngs flow (to first order) as 

dgk/dl = {d- Xk)gh - X] (^ijk9i9j + ■■■ 

ij 

where Cijk are the OPE coefiicients 

Before the randomness is introduced, the operators Sc^d • • • with q distinct rephca labels are 
degenerate. We need to calculate the OPE coefficients of the disorder operator O = '^a^h^o-^b 
within this subspace. In SU (2) language let us denote the presence of an energy operator in replica 
i by Wi) and its absence by | J.i). We can write a general operator by the action of a series of SU{2) 
raising operators on a 'vacuum' |J,i|2 • • • in): 

ea£b ■ • • < > T+T^ . . . |iii2 ■■■In) 

with 

1 " 1 
g=-^(l + 2rf) = -n + ^rf. 

i=l i 

In the subspace corresponding to the energy sector, the OPE of one of these operators with the 
disorder operator O is equivalent to the action of the matrix 

on the corresponding state, with [-Mj^^rf] = as we would expect. (The extra factor of two 
arises from the two different ways of performing the operator contraction). If we define the total 
spin vector S = X)iLi then we can write 



i 

= 2{S^ - {S'f + - q) 
Using 5 = 5^ + § we can thus write 

M = 2(s{S + l)-{S^f-'^). 

However, we now have the problem of interpreting the total spin angular momentum S. One 
solution is to ask which values of S are consistent with the known value of — clearly, the 
possible values are 

_ n n n n 
2' 2 ~ ' 2 ~ ' ' 2 

We can support this by comparing the irreducible representations of SU (2) with the irreducible 
representations of Sn considered by Ludwig. In the Sn case, we have to consider the representations 
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given by all Young tableaux of n boxes and at most two rows, with up to q boxes in the second 
row. E.g., for g — 2, n ~ 6 the possible irreps are pO| ] 



e 


e 













These irreps can be placed in one-to-one correspondence with the representations formed by 
the tensor product of n spin-half representations of SU{2) with precisely q of them spin-up. Of 
these latter tableaux, a tableau with (n — r) boxes in the first row and r boxes in the second row 
has a dimensionality of D = n — 2r -|- 1 corresponding to a total spin of 5' = ^^j^ = ^ — r, as was 
proposed by the consistency condition. Hence, in general the q*'' moment of the energy operator 
will have a series of irreps 3q^r with scaling dimensions 

where r G {0, 1, . . . ,q}, y is the RG eigenvalue of the coupling to disorder and is the scaling 
dimension of the energy operator at the pure fixed point. 

The factor of j arises from the position of the disordered fixed point, set (to one loop) by 

dA/dl = aA + 4(n - 2)A'^ + 0{A^). 
After taking the n — > limit we obtain 

X,AO) = 9 - f {rir -l)-q'+q)+ 0{y^). 

The most relevant scaling dimension belongs to the most antisymmetric irrep 3q^q and has a linear 
dependence on q. Note that the irreps 3q_o s-^d ^g.i ^re always degenerate at n = 0. This must 
happen because for n ^ they have degeneracies 1 and (n — 1) respectively, while as n ^ all 
operators must arrange themselves to have degeneracies oc n, so that the torus partition function 
is 1 -|- 0{n). In fact, this collision of the scaling dimensions at n = will give rise to a logarithmic 
operator 

The operators ei . . .£q span a vector space of dimension (p. In general, a Young tableaux of 
shape [n — r, r] corresponds to a vector space of dimensionality ^ shown 

that the sets of permitted irreps (r G [0,(7]) exhaust this vector space. The dimensionality of a 
given irrep does not depend on the number of energy operators in the tableau, and is preserved 
under the addition of extra e-operators via — irreps related via the action of are joined in 
the plot below (figure [|) : 

The two-point function of the most antisymmetric irrep for a given g, corresponds to the 
quenched average of the q^^ moment of the connected correlation function, 

Gl = (ei - £2) ■ • • (£29-1 - £2g)(0)(£l - £2) . . . (£29-1 - £2q){r) ■ 

This can be seen by noting that the state corresponding to the operator (£1 — £2) . . . (£29-1 — £2q) 
is annihilated by the step-down operator S~ = '^^ '■ 

(E (^1+ - r+) . . . (r^Vi - r2+) Uii2 ... in) = 0. 

All T~ with i > 2q give zero when acting on everything to their right, so we only need consider 
those T~ with i < 2q. Take the terms of the step-down operator in pairs: the first pair (ji + t^) 
will act on everything to the right to give 

(rrr+ - t,t+){t+ - T+) . . . (r+_i - r+) [UU ■ ■ ■ in) ^ 0. 
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Figure 1: Scaling dimensions of q moments of the energy operator to one loop at y = 0.1 



The action of all other pairs similarly vanishes, so the state corresponding to the connected cor- 
relation function is annihilated by and is therefore also the state corresponding to the most 
antisymmetric irrep. 



3.2 Mixed Sector: Exclusion argument 

We now wish to consider the behaviour of cross-moments, i.e., investigate the renormalisation of the 
operators ai . . . UpSi . . .Sq. Fortunately, given the behaviour of the ei . . .Sq operators considered 
above this is not a particularly difficult extension. To one loop, we again consider the effect of 
the disorder operator J^a^^b ^a£b on the operator . . . ej^ . . . Sj^ , and look for operators of the 
same form to be produced by all possible contractions. These contractions are given by the fusion 
rules of our theory at the pure fixed point: 

£o • £6 Sa,b 

Oa-Ob ~ ^a,b + + '^(Q - 2)) Ea 

Contractions are only possible within a given replica. Given the fusion rules, at this order 
there are only two possibilities: 

1. a e {j}, 6 ^ {i\ U {j}: one of the energy operators in the disorder operator contracts against 

another energy operator {ee ^1), the other lies in an 'empty' replica. This gives the same 
behaviour as in the pure energy sector, but with a shift in the effective number of replica 
indices: n — > n — p 

2. a, 6 € {i}, a^h: both energy operators in the disorder operator contract with spin operators, 
giving rise to further spin operators. As considered by Ludwig, this gives rise to a term 

\P{P - 1) 

Any other attempted contractions will evaluate to zero when the the trace over the fields is taken. 
Substituting these changes into our expression for the energy sector, after taking the n ^ limit 
we find 



where as before r G [0,(7] and is the scahng dimensions of the spin operator at the pure 
fixed point. The spin operators exclude the energy operators from their rephcas, giving rise to an 
effective shift in the number of replicas n — > n — p. 

The scaling dimensions of the most antisymmetric irreps r = q are unaffected by this ex- 
clusion process, as they have no n-dependence. This means that, to one loop, there will be no 



extra structure in the quenched average {a{Q)<T{R)) {e{0)e{R))^. In section 3.4, however, we shall 
demonstrate that this quantity does have extra structure for the case of the random-bond Ising 
model in 2 -|- e dimensions. 

By analogy with the SU(2) argument, it is possible to consider the mixed sector by writing 
our operators {e, a, 1} as a fundamental representation of SU (3). This analysis leads to the same 
results, as the {l,e} SU{2) subgroup effectively decouples from the spin operator. 

As one check on these results, standard arguments in probability theory (see, e.g., [T^ ) require 
that —Xp q q should be convex for all increases in p or q, both individually and jointly. Considering 
the second derivatives of this quantity w.r.t. p and q, this requirement can be seen to be trivially 
satisfied for y > 0. 

3.3 Replica Structure 

Given that non-trivial structure in the spin-energy sector only appears to sub-leading order, it is 
useful to be able to write two-point functions of irreducible representations of the replica group 
in terms of numerically available quantities. We have already noted that the most relevant, most 
antisymmetric irreps correspond to moments of the connected part of the two-point correlation 



function (section 3_^) and would like to be able to say something similar for a general irrep. In 
Appendix A we detail the combinatorics needed to decompose a two-point function of irreps into 
a linear combination of various other two-point functions: here we shall merely note a few results. 
We find explicitly that the most antisymmetric irrep at a given 3q^q has a two-point function 

{Jq^q{0)3q.q{R)) ^ {{Ei - £2) ■ • • (£2^-1 ^ £2q){0) (EI " £2) • • • (£2?-! " £29)(^))rop 

In the q = I case, both the possible irreps Ui.o and correspond to the connected correlator 



(e(O)e(i?)} — (£(0)) (£(i?)). This is in fact a general feature, since the lowest two irreps for arbitrary 
q are always degenerate after taking the n — > limit. 

As an example of a correlation function which has a non-trivial behaviour on addition of p 
spin- spin correlators, the irrep J 2,1 is a non-leading irrep which can be written as 



(J24(0)32,l(i?)) - (£(0)£(i?))(£(0)£(i?)) - 4(£(0)£(i?))(£(0))(£(i?)) + 3(£(0)) (£(0)) (£(i?)) (£(i?)) . 

If we take a factor of {a{Q)a{R)y inside each quenched average on the RHS, the resulting scaling 
dimension will be 

Xp,2,i - pxa + 2x, - |(-8 - 2p + i(p2 -p)) + 0{y^) 

^P,o,o + ^0,24 = + 2x, - |(-8 + i(p2 - p)) + 0{y^) 

Thus this correlation function directly exhibits the cross-structure between the two-point function 
of the spin operator and the full two-point function of the energy operator. 

3.4 Ising model in 2 -|- e dimensions 

As an alternative to perturbing the pure Ising model by changing the number of possible values 
for the spin, we can consider the Ising model in 2 -|- e dimensions. The specific heat exponent a is 
zero ai d — 2 and small for d = 3 (< 0.1), so we shall assume that a = 0{e) with e ^ 1. The bond 
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disorder is relevant for all e > 0, and we assume that the system will flow to a nearby disordered 
fixed point. The calculation of scaling dimensions at this fixed point is very similar to that for the 
disordered Potts model — the only difference is that the OPE rules are altered to 



5afi{l + ce) 

,5,,b(l+ ((l/2) + 0(e))e) 



The OPE coefficient c is zero for d = 2 due to the self-duality of the Ising model in two dimensions 
1^, chapter 8], but non-zero for d > 2. Repeating the analysis of section 3.2 we can obtain two 
more types of contraction between the disorder operator X^a^^b ^a.£b and aPe'^ : 

1. Both energy operators in the disorder operator contract with energy operators, giving rise 
to two further energy operators in the same replicas. This gives a term c^q{q — 1). 

2. The disorder operator contracts with one energy operator and one spin operator, giving 
rise to one operator of each type with unchanged replica indices. This produces a term 
2 • (1/2) ■ c ■ pq — cpq. 

The non-zero c will shift the position of the fixed points, via the RG equation for the disorder 



dA/dl = a A 
so we obtain scaling dimensions 



(4(n - 2) + 2c2) + 0(A^ eA^ 



^p,q,r = PXcr + qXe + 



y 



p 



+ ^{P -P)+ c\q' -q) + cpq\+ 0{y') 



The cpq term indicates that, with a non-zero c, we now have cross-structure between the two- 
poi nt function of the spin o pe rator and the connected tw o-point function of the energy operator, 
i.e. {a{0)a{R)r {e{0)e{R))l + {a{0)a{R)f ■ {s{0)s{R))l 



4 Conclusions 

To one loop, the interaction between the spin and energy operators is an 'exclusion effect': the 
spin operators block replicas, shifting the n-dependence of the scaling dimensions of the different 
irreps in the energy sector. For the Potts model this does not have any effect on the leading 
behaviour, as the most relevant irrep irrep r = q has a scaling dimension with no n-dependence. 
The sub-leading terms are affected, however, and this may be picked up in numerical studies. By 
writing the irreps in terms of correlation functions we have shown how linear combinations of 
correlators could be used to demonstrate the existence of an underlying joint distribution function 
for the quenched average of spin/spin and energy/energy moments. In the Ising model in (2 + e) 
dimensions the leading behaviour of mixed moments is affected non-trivially. Also, three- or 
higher- point functions should also exhibit this cross-sector behaviour, presumably with selection 
rules coming from the group structure of the energy sector. 

There remains the question of what happens when these calculations are continued to two loops. 
A priori, there seems no reason to expect that the scaling dimensions of the most antisymmetric 
irreps will continue to be protected against corrections from cross-correlations between the spin 
and energy two-point functions. We will content ourselves here with noting that a coulomb-gas 
calculation along the lines of would give two-loop corrections to the scaling dimensions of the 
operator ae proportional to 

J d'y {e{0)e{y)e{l)e{oo)) {e{0)e{y)<j{l)<j{^)) 

which we would expect to be non-zero, although a detailed calculation would of course be necessary 
to confirm this. 
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A Decomposition of {'3q^r{0)'^q,r{R)) into quenched averaged 
energy correlators 

Start with the most antisymmetric irrep for q energy operators, "^q^q, corresponding to the Young 
tableau 



El 




Sq 










n — 2q 



The correlator {3q^q{0)3q^q{R)) is proportional to 

n — q 

(ei — e„)(£2 — £n-l) ■ ■ - {Sq— e„-q+l)(0) ^ {Sa' — " ^n-l) ■ • ■ (Cg' " Sn-q+l){R)- 

a'^b'^...^q' = l 

We shall start by expanding this into 2^'^ monomials: 

• Choose N operators from the LHS with index G [n — q + l,n], < N < q. There are 
J2n=o (n) '^ays of doing this. 

• Choose M matching operators on the RHS, < M < A^: J2m=o (m) ways. 

• Choose O other operators with index G [n—q+1, n] on the RHS, < O < q-N: Y^o^a {"'J^) 
ways. 

Note that X)w=o (n) Sm=o (m) X]o=^ (^o^) ~ 2^'^' have accounted for all the terms. 

For a typical monomial 

El . . . Eq-NSq ■ ■ ■ £q-N+l £a' ■ ■ • £((;-A/-0)' • ■ ■ ^q-M+l^q-N ■ ■ ■ Eg-JV-O+li 

the number of explicit free indices ej^e'^, ... on the RHS that are paired, P, runs over the range 
< P < Mm{q — N,q — M — O). For a given P we have the following factors: 

• A sign of {-l)N+M+0 

• The number of ways of choosing the pairings, ("^p^) ('^^^^*^)P! 

• {q~ M ~0 ~ P) explicit, unpaired free indices ranging over [n — q) — (q — N) replicas, giving 
a factor of 

{n + N -2q)\ 
{n + N -Sq + M + + Py. 

• (M + O) implicit free indices on the RHS (which must still be summed over) ranging over 
{n — q) ~ {q ~ M — O) replicas, a factor of 

{n-2q + M + 0y. 
{n-2qy. 

A term with P pairings will produce a correlator 

Cp+M, q-p-M - {e{0))<i-P-^' {e{Q)e{R))P+^i (e(i?)}?-^-^^ 
Altogether, our decomposition into correlators becomes 

^=0 ^ ' M=0 ^ ' 0=0 ^ ' P=0 \ / \ / 

,,jv+M+o (n + jV-2g)! (n-2g + A/ + 0)! ^ 

^'^^ („ + iV-3o + M + + P)! ^7^2^M ^P^M,q^P-M 
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For various q, this formula gives 



q=l 
q = 2 
q = 3 

q = A 



n{Cifi - Co,i) 
(n - l)(n - 2)(C2,o - 2Ci,i + C7o,2) 

(n - 2)(n - 3)(n - 4)(C3,o - 3^2,1 + 3Ci,2 - C„,z) 
(n - 3)(n - 4)(n - 5)(n - 6)(C4,o - 4^3,1 + 6^2,2 - 4Ci,3 + Co,4) 



g = 5 (n - 4)(n - 5)(n - 6)(n - 7)(n - 8)(C5,o - 5C^4,i + lOCa.z - lOCa.a + 5Ci,4 - Cq.s) 

i.e., the most antisymmetric irreps correspond to moments of the connected two-point correlation 
function Gt- 

The above argument can easily be extended to irreps which are not the most antisymmetric. 
If we consider now the irrep 3q_r corresponding to Young tableau 



e 




e 


e 




£ 










• V V ' 

q—r n—q—r 



the following changes occur: 

• The {q — M — O — P) explicitly unpaired indices now range over (n — r) — (g — N) replicas, 
giving a factor of 

{n + N -q-r)\ 



{n + N -2q-r + M + + P)\ 

• The (M + O) implicit free indices range over {n — r) — {q — M — O) replicas, a factor of 

{n-q-r + M + Oy. 
{n — q — r)\ 

Incorporating these changes, we obtain (for instance) 

g=l,r = Ci,o + (n - l)Co,i 

g = 2, r ^ 2C2,o + 4(n - 2)Ci,i + (n - 2)(n - 3)Co,2 

5 = 2, r = 1 n(C2,o + (n - 4)Ci,i + (3 - n)Co,2) 
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